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Abstract 

A set of N three component unit scalar fields in (3 + 1) Minkowski 
space-time is investigated. The highly nonlinear coupling between 
them is chosen to omit the scaling instabilities. The multi-soliton 
static configurations with arbitrary Hopf numbers are found. More- 
over, the generalized version of the Vakulenko-Kapitansky inequality is 
obtained. The possibility of attractive as well as repulsive interaction 
between hopfions is shown. A noninteracting limit is also discussed. 
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1 Introduction 



Topological defects are one of the main ideas of the temporary theoretical 
physics. They play crucial role in many models - from cosmological strings 
and textures pQ, [2] to magnetic monopoles and vortices in Quantum Chro- 
modynamics jl]- On the other hand, topological defects are observed in 
various condense matter experiments - see for instance liquid crystals and 
3 He [H] or 4 He quantum liquids [Hj. Among quite well understood topolog- 
ical defects like domain walls, strings/vortices and monopoles toroidal-like 
configurations, that is solitons with nonzero value of the Hopf index, seem 
to be rather mysterious objects. 

Since Faddeev, Niemi [7], jH] and Cho [5] have proposed their famous effec- 
tive model for the low energy gluodynamics, where glueballs (particles which 
consist of only the gauge fields) are knot-like objects made of self-linking flux- 
tubes, toroidal solitons, their physical as well as mathematical aspects have 
been considered by many theoreticians l . Unfortunately, due to the rather 
complicated toroidal symmetry and nonlinearity of equations of motion no 
analytical results, in recently investigated QCD relevant models [12], [T3j . 
|14j . have been found. Almost all results have been obtained by means of 
numerical methods [TH], [IE!, C3> P^l - I n or der to deal with exact toroidal 
solutions one is forced to consider various toy models. 

All such theories, providing analytical description of hopfions, are based on 
the quite old idea presented by Deser et al. [19,. They have investigated a 
classical highly nonlinear field theory with a form of the Lagrangian chosen 
to circumvent Derrick's theorem. In fact a soliton solution with unit Hopf 
index [20] and then configurations with arbitrary value of the Hopf invariant 
have been found [21]. Moreover, nontrivial topological solutions have been 
obtained in more complicated model [22], for example with explicitly broken 
0(3) symmetry [23] . Using the exact soliton solutions one has checked that 
the Vakulenko-Kapitanski [21] inequality is fulfilled in all these models. 
It should be mentioned that soliton solutions with non-trivial Hopf index 
can be achieved also in another famous model i.e. Skyrme model in (3 + 1) 
dimensions [23], j2E]- Moreover, there is also possibility to analyze hop- 
fions in less dimensional space-time. It is based on the observation that the 
Faddeev-Niemi knots can be understood as twisted magnetic vortex rings 
built of so-called new baby skyrmions found in (2 + 1) anisotropic Skyrme 
model [23- 

In the present paper we would like to investigate toroidal soliton solutions in 
1 Knotted solitons appear also in condensed matter |10| and possibly in astrophysics 
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a model which consists of many unit vector fields i.e. allows for many topo- 
logical charges. The appearance of more than one unit vector field has been 
observed in recently derived non-Abelian generalization of the color dielectric 
model ^3]. Indeed, there are two unit fields in this model. The main aim, 
besides finding exact form of the topological solutions and their energies and 
Hopf charges, is to analyze Vakulenko-Kapitanski formula in case of more 
than one Hopf index. 



2 The model 

Let us consider a set of N unit three component scalar fields Hi, i = 1...N, 
H\ = 1, living in the (3 + 1) Minkowski space-time. The Lorentz and 0(3)^ 
invariant Lagrangian density is chosen in the following form 

N 



c= n (fa ■ x d ^i)f) at 



i=l 



which is a generalized version of the Aratyn-Fereira-Zimerman model with 
N unit fields. It can be rewritten as C = Y\f=i{H^} 2 ) a \ where = Hi ■ 
(d^Hi x d v Hi) are antisymmetric field tensors. In order to omit the Derick 
arguments for nonexistence of static soliton solutions one has to assume the 
condition 

N o 

E«i = 7- ( 2 ) 

In fact, now the total energy becomes invariant under the scale transforma- 
tions. 

Now, after taking advantage of the stereographic projection, we can express 
each of the unit fields in terms of two complex functions u and u*: 

Hi = — -, — 7^{ui + u*, -i(ui - u*), \ui\ 2 - 1). (3) 

1 + \Ui\ 2 

Then the Lagrangian density reads 

TV 1 

^ = 8 3/4 n (1+k|2)4Ql (^xr% (4) 



where we have introduced a set of the objects: 

Kf = (d^Ui^m - {duUicTu^u*, (5) 
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where i — 1...N. It is straightforward to check that they fulfill the following 
conditions 

Kfd^m = and Im(Kfd»u*) = 0. (6) 
The equations of motion take the form 



N 

n 



:i + \ui\ 2 )^ 



N 



n 



; (1 + kP) 4Q < 



1 + M 2 ) 4 ^' 
1 



(i + 



(7) 



After some calculation one can rewrite them as follow 



0„ 



n 



or in the more compact form 



1 + H 2 ) 4a » I (1 + \uj\ 2 Y a r 



ax (j> = o, 



where 



n 



;i + | M .|2)4a, I (l + | Uj .| 



(9) 



(10) 



It is easy to show that these quantities also fulfill the previously found con 
ditions (JUJ). Because of that we can generalize the procedure defined in 
and construct the infinity families of the conserved currents 



J. 







dui du* ' 



(11) 



where Gi, i = 1...N are any functions of Uj and u*, j = 1..N. Due to that 
presented 0(3)^ invariant model is integrable in the sense that N infinite 
families of the conserved currents can be found (2§j- 
Let us proceed further and prove that, as in case of soliton models in (1 + 1) 
and (2 + 1) dimension, the existence of infinite number of the conserved 
currents leads to topological soliton solutions. We begin with introducing of 
the toroidal variables 

x 



q 



sinh rj cos ( 
sinh rj sin < 
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a ■ c 
-sinf, 

q 



(12) 



where q = cosh 77 — cos £ and the constant a > sets the scale of the coordi- 
nates. However, one should keep in mind that the model is scale invariant. 
Due to that, there is no scale on the solution level. All solitons, irrespective 
of size, have identical energy. 

In addition we assume a natural generalization of the Aratyn-Feriera-Zimer- 
man Ansatz 



(13) 



where rrii,ni, i = 1...N are integral parameters. 

Inserting the Ansatz into the field equations we derive at the following 
static equations 



d v In 



N 

n 



{hfl 



l\2ai 



(./:,/; 



2\ 4at 



;i + f 



-d v In 



N 



sinh 77 Y[ ntf H — : 



i=i 



sinh ?7 / 



These equations can be integrate and we obtain that 



(14) 



N 

n 



N 



1 + /, 



1 + / 



^3 

sinh 77 



2 

r?.f 



i=i 



sinh 77, 



(15) 



where the integration constants fcj, j = 1...N have been introduced. 

One can show that (|T5j) can be reduced to the set of N decoupled first order 

differential equations 



1 



N 



4Ql m 2 + 



7?T 



sinh 77 / 



sinh 2 77 ' 



(16) 



with the following general solutions 
1 1 



f 




4q < m 2 + 



??7 



-2a> 



sinh 77 



sinh 77 



(17) 



Here /j, j = 1...AT are integration constants. 

Now, we will analyze the total energy corresponding to our model (|TJ): 



E j d 3 xT 00 = 8 3 1 j r/ :, r J] jj-' ~ 



i=l 



+ \Ui 



2^4^"' 
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where the stereographic projection © has been taken into account. Inserting 
(0) and the Ansatz (JT3j) into (JTSj) we derive that 

£m ,„ = ^ . (n ('-■ + AD (19) 

Quite surprisingly, using the equations (fT7)|) we are able to remove the un- 
known functions fj from this formula. Then the total energy integral takes 
the form 

, N f°o dn N ( n 2 \~ 2a * 

E m , n = ( 27 r) 2 8-2 3 / 4 n^7 ^IlK + ^r , (20) 

J o smh T] " \ smh r\ ) 

which can be evaluated, at least by means of some numerical methods, for 
all values of the parameters rrii, Hi, ccj. 

One can notice that behavior of the integral (and difficulty of the calculation) 
strongly depends on a particular value of the following ration 

S - * < 21 > 

In the next section the simplest case, when all these numbers are equal will 
be investigated. 

3 Case q 2 = const. 

Let us focus on the case when 

Tl 2 

_L = q 2 = const., (22) 
mf 

for all i = 1...N. Then the equations of motion (JTfij) can be simplified to the 
following expressions 

(i+% 2 = kj n ^ ( m * 2 + ^n^) s^v (23) 

Moreover, these equations can be integrated out and finally we obtain that 



1 + // kjl-q 



1 1 ^tt,4q -4«-\ cosh 77 , . ^ 

' n C'mi 4Ql J ~~r + /j, (24) 

g 2 + sinh 2 77 ) 2 



2 



vi=l 
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where j = 1...N. Here we have taken q 2 < 1. In order to find the value of 
the integration constants we need to specify the asymptotical conditions for 
the functions fj. We choice them in the form which admits the nontrivial 
topological structure |2I] 

(0,0,1) i.e. / -> oo as rj -> (25) 

and 

n -> (0, 0, -1) i.e. / -> as r? -> oo. (26) 

Then we find that 

nf =1 tf* i i - g 2 



kj " 2Ilf =1 m- 4Q ' |g| - 1 



\q\- (27) 



This is equivalent to the fact that all integration constants have to take the 
same value 

kj = k = const. (28) 
In addition, one can obtain that 

l 3 = l = --^—r- (29) 
\q\ — 1 

Finally, the solution is given by the following formula 

1 1 



1 + ff l~\q\ 



cosh i] 

1 - \q\ 



q 2 + sinh 7] 



(30) 



One can recognize here the set of N Aratyn-Fereira-Zimerman toroidal solu- 
tions. 

Let us now show that the solution possesses toroidal symmetry. The compo- 
nents of the unit vector fields read 



n h = i r V 2 cos ( m i£ + n i4>), (si) 



n (>) = 1 + f? Sill ( m ^ + n ^)' ( 32 ) 

3 ~ 1 + fi 
U (i) = 1 , A ' (^S) 



where /, is given by (}3T)j) i.e. depends only on the radial-like coordinate rj. It 



is clearly visible that surfaces of constant nf^ are torii. 
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After inserting obtained here solution into the total energy integral one finds 
that 



N i i— 

E m , n = (2vr) 2 4 • 2 1 ' 4 II m^ijl + \q\y/\q\. (34) 
i=i 

In the limit m ai = m, i = 1...N the formula founded in j2l] can be repro- 
duced. 

In order to prove that our solutions possess non-trivial topology one has to 
calculate value of the pertinent topological invariant. In our case static field 
configurations, which can be viewed as maps from S 3 into S 2 , are classified 
by so-called Hopf index Q H = n 3 (S 2 ). One should notice that because of the 
fact that we deal with N fields thus N different topological Hopf charges 
have to be introduced. 
We define two additional functions 



n 2 



g\ = cosh 77 — \l + sinh 2 r] (35) 

and 



m 2 



2 / 171 2 

g 2 — ylH T sinh // — cosh?]. (36) 

V n z 

Then, for every field we introduce four objects 



2 



and 



'9i+92. 




= ~r= >< ftnt) • (3? 



They are related with the i-th unit vector field via the standard relations 
Hi = Z^ l aZ l , where a are Pauli matrices and 



^2 



= (Z\* , Z l 2 *) (39) 



and 

Z\ = + i<f>®, Z\ = + i$f. (40) 

Using these objects one can construct N Abelian vector fields which form 
the previously introduced antisymmetric field tensors = dkA® — fyA®. 
In fact one finds that 

Af = % -(Z i *d k Z i -d k Z i *Z i ). (41) 
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Then, using the standard expression for the Hopf index 

Q$ = ±iJ<Px&>-m, (42) 
we arrive at the following result 

= -riirrii. (43) 

In other words our soliton solutions are classified by set of N Hopf charges. 
Every kind of the unit vector field carries its own topological charge which is 
unique determined by the parameters raj, rij. Now, one could ask about gen- 
eralization of Vakulenko-Kapitansky inequality for N vector fields. Knowing 
that yj\ + \q\y/q > V2\q\ 3 ^ we obtain the generalized relation 

N 

£>(27r) 2 4-2in^)' (44) 

i=l 

where J2iLi a i — f • If & U topological numbers are equal Qrq = Q, i = 1..N 
we reproduce the standard Vakulenko-Kapitansky formula. 



Qh 



(9\2 



2\2 



(($?: 



(i)\2\2 



4 Conclusions 

In the present work a generalized Aratyn-Ferreira-Zimerman model has been 
analyzed. This model consists of N unit vector fields and possesses N inde- 
pendent 0(3) symmetries. We have proved that this model is integrable in 
the sense that infinite number of the conserved currents appears. Precisely 
speaking, one can find a set of N infinite families of the currents. 
We have also shown that, as in the case of the soliton theories in (2 + 1) 
dimensions, this property leads to the existence of non-trivial topological 
field configurations. Our soliton solutions are described by shape functions 
depending on rj coordinate and angular terms characterized by set of two 
integer parameters raj and rij. The total energy has been also found. In gen- 
eral, it is given by the integral (|20j) . The exact form of the shape functions 
has been obtained in the limit of the constant ratio of the Ansatz parameters 
nil ml = q 2 = const. Then the shape functions take form of the well- know 
Aratyn-Ferreira-Zimerman solution. However, obtained solutions are much 
more general due to the angular part which for each of the unit vector field 
can be different. Because of that the topological behavior of the solutions is 
changed - now we have N independent topological Hopf indexes. In other 
words we have obtained multi-soliton configurations classified by iV topolog- 
ical numbers Q^ 1 , where each of the vector field carries different topological 
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charge. These configurations consist of N toroidal solitons of the different 
types. 

It should be stressed that all properties of the solutions (topological charge, 
total energy), at least in the qf = const, determined only by value 

of the parameters m^nj. For example, behavior of the shape functions fi 
depends only on the ratio between them whereas number of the fields does 
not affect the shape functions. Regardless of the particular number of the 
fields N, all functions fi possess identical form. On the other hand, number 
of the vector fields defines the topological content of the model. Similarly, 
particular values of the parameters on of the discussed model also do not 
inflect on the shape of solitons as well as on their topological charges. 
In our opinion, the most interesting result is the generalized Vakulenko-Kapi- 
tanski inequality. As we see this formula can be divided into two parts. The 
first part contains only numbers and seems to be universal i.e. Lagrangian 
independent. In fact, it is identical as in the simple one vector field case. The 
second part, which includes the Hopf indexes, appears to be strongly sensi- 
tive to the form of the model. To show it more precisely one can consider a 
simple sum of N standard Aratyn-Ferriera-Zimerman Lagrangians 

N 

£ = £(#if) 3/4 

i=l 

i.e. a model with the noninteracting unit vector fields. One can immediately 
show that the pertinent Vakuleko-Kapitanski formula reads 

TV 

E > (2vr) 2 4-22^|Q w |i. 
i=i 

We see that the second part of the inequality reflects a form of the La- 
grangian. However, a common feature is shared by both inequalities - there 
is an universal scaling property. Namely, if we transform all topological 
charges — > then the total energy will scale as E — > X*E. 
Let us now analyze the physical meaning of the parameters Oj. In spite of 
the fact that the form of the solution, at least in the qf = const, regime, 
does not depend on their values, they play crucial role in the interaction of 
the toroidal solitons. Indeed, interaction between hopfions in the proposed 
model is governed by the values of a, and it can be attractive or repulsive. 
Quite interesting hopfions can even do not interact at all. To be more pre- 
cisely we will consider two hopfions with the following topological charges 
Q a = (Q(i) ? .., QW), a = 1, 2. It follows from the formula that if a { < 1 
then i-th hopfions built of the z-th vector field should attract each other 
whereas for a, > 1 this interaction will be repulsive. In the very special 



10 



case where Oj = 1, they do not interact at all. It resembles the Bogomolny 
limit in the standard soliton systems. As there is only one constrain on the 
parameters ai (j2j), all three situations can be realized. Thus, in general, one 
can deal with a model where every possible type of interaction occurs. 
As we have mentioned it above, the analytical solutions have been obtained 
only in the simplest situation i. e. for qf = const.. It constrains also the 
validity of the generalized Vakulenko-Kapitanski inequality. Thus one should 
investigate more complicated cases with qf ^ q? as well. We would like to 
address this issue in our next paper. 

This work is partially supported by Foundation for Polish Science and ESF 
"COSLAB" programme. 
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